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X-ray absorption spectroscopy is a well established tool for obtaining information about orbital 
and spin degrees of freedom in transition metal- and rare earth-compounds. For this purpose usually 
the dipole transitions of the L- (2p to 3d) and M- (3d to 4f) edges are employed, whereas higher 
order transitions such as quadrupolar Is to 3d in the K-edge are rarely studied in that respect. This 
is due to the fact that usually such quadrupolar transitions are overshadowed by dipole allowed Is 
to 4p transitions and, hence, are visible only as minor features in the pre-edge region. Nonetheless, 
these features carry a lot of valuable information, similar to the dipole L-edge transition, which is 
not accessible in experiments under pressure due to the absorption of the diamond anvil pressure- 
cell. We recently performed a theoretical and experimental analysis of such a situation for the metal 
insulator transition of (V(i_aj)Crcc)203. Since the importance of the orbital degrees of freedom 
in this transition is widely accepted, a thorough understanding of quadrupole transitions of the 
vanadium K-pre-edge provides crucial information about the underlying physics. Moreover, the 
lack of inversion symetry at the vanadium site leads to onsite mixing of vanadium 3d- and 4p- states 
and related quantum mechanical interferences between dipole and quadrupole transitions. Here we 
present a theoretical analysis of experimental high resolution x-ray absorption spectroscopy at the 
V pre-K edge measured in partial fluorescence yield mode for single crystals. We carried out density 
functional as well as configuration interaction calculations in order to capture effects coming from 
both, itinerant and atomic limits. 

PACS numbers: 71.27.+a, 71.10.Fd, 71.30.+h 



I. INTRODUCTION 

The process of X-ray absorption spectroscopy (XAS) 
is the excitation of a core-shell electron under specific 
selection rules. Experimental techniques evolved, to- 
gether with the theoretical understanding and simula- 
tion capabilities over the last 25 years to maturity and 
have hitherto proven to be a valuable tool for probing 
the groundstate of transition metal (TM) and rare earth 
(RE) compounds. Advanced techniques make use of the 
polarization dependence of the absorption. Measuring 
the dependence of the absorption for circularly polar- 
ized light, called circular dichroism (CD), together with 
powerful sumrules provides a spectroscopic way of mea- 
suring the magnetic susceptibility^'^. The dependence of 
the absorption for linearly polarized light, called linear 
dichroism (LD), on the other hand, has been success- 
fully used to detect the orbital occupation in various sys- 
tems such as high Tc cuprates^, Ce-based heavy fermion 
compounds^, and also oxides of vanadium like VO2 and 
V2O3 (Refs. iandi). 

For the transition metals the most informative excita- 
tions are obviously those involving 3d orbit als. In this re- 
spect, most of experimental and theoretical analysis have 



been done for the V2p ^VSd transitions, the so called L- 
edges (at energies around ~ 350 eV to ~ 950 eV - i.e. soft 
x-ray range). These transitions are dipole allowed and 
therefore have large cross-sections. Nonetheless, com- 
plementary valuable information can be extracted also 
from other core electron excitations of higher order like 
Vis ^V3(i quadrupole transitions located in the pre- 
edge region of the so called K-edge (at energies ranging 
from ~ 5 keV to ^ 10 keV for transition metals - i.e. 
hard x-ray range). 

One of the remarkable advantages of XAS is the possi- 
bility of an element selective measurement^, since the ab- 
sorption edges have distinct energies due to the element- 
specific binding energies of the core shell electrons. Im- 
portant for the application of X-ray absorption as a spec- 
troscopic technique is the development of models to simu- 
late accurately the observed spectra: a quantitative anal- 
ysis allows for the determination of various (near) ground 
state properties, including the valence, spin and orbital 
state of the atoms under investigation^. These calcula- 
tions are based on cluster models and have been success- 
fully extended^^^ as to include the influence of the core- 
hole using the full atomic multiplet theory. The scheme 
relies on the localized nature of the exciton and consists 
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basically in diagonalizing a configuration interaction (CI) 
Hamiltonian. The calculations are parameter-based and, 
hence, not ab initio. Mostly, the parameters are fitted to 
experimental data in order to extract information. The 
"cluster calculations", as they are often called, became 
quite popular (especially for the L- and M-edges) due to 
their - sometimes outstanding - ability to reproduce ex- 
perimental data. Ab initio methods, on the other hand, 
such as density functional theory in the local density ap- 
proximation (LDA) often reproduce well the continuum 
part of the spectra or the main edge region of K -edges. 
However, the density functional methods fail when exci- 
tonic features become important. Hence, in the broadest 
sense the two methods are either tailored for the limit of 
strong (atomic) electronic correlations (CI) or the itiner- 
ant limit (LDA). 

With this in mind the TM K-pre-edge region^ is pecu- 
liar since it inherits a "meeting point" of the atomic- and 
itinerant limit which, by its own right, makes such a chal- 
lenging analysis appealing. Yet, this is not the only mo- 
tivation to obtain a better understanding of higher order 
transitions in pre-edge regions of TM K-edges. In fact, 
the main theoretical interest for the K-pre-edge transi- 
tions originates from the possibility of analyzing XAS 
experiments performed under external pressure: while 
the use of external pressure provides one of the clean- 
est ways to change the physical parameters of a com- 
pound, in such experiments the dipole allowed L-edge 
transition energies overlap with the energy scale of the 
absorbtion of the diamond pressure cell, so that the K- 
edge in the hard x-ray regime represents a technique 
which can provide information inaccessible with other 
methods^^. This problem is a common one when phase 
diagrams/transitions are explored by means of applying 
pressure. And, naturally, it becomes a severe obstruction 
when orbital or spin degrees of freedom play a major role 
in the studied compound. A well known example of a 
pressure induced transition is the metal-insulator tran- 
sition (MIT) in chromium doped vanadium sesquioxide 
(V(i_2;)Craj)203 J^ii^ First considered to be a textbook 
example for a simple single-band Mott insulator, it was 
later realized that the metal insulator transition upon 
doping pure V2O3 (metal to insulator) with chromium 
or applying pressure to the doped sample (insulator to 
metal) cannot be understood without incorporating cru- 
cial orbital physics^ii^"— . While this issue of orbital oc- 
cupations was studied for the doping-induced transition 
a while ago employing the vanadium L-edge^, a thorough 
understanding of the pre-edge region in the vanadium K- 
edge is still missing, hampering the analysis of the pres- 
sure induced transition. In a recent work2^ we presented 
first experimental data for the vanadium K-edge as a 
function of temperature, doping, and prssure in a pow- 
der sample (for an extended experimental report see also 
Ref . [23I ) . With a theoretical analysis by means of LD A+ 
dynamical mean field theory (DMFT) combined with CI 
cluster calculations we pointed out a clear experimental 
evidence showing the difference between the metal insu- 
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FIG. 1. (color online) Conventional (left hand side) and primi- 
tive (right hand side) unit cell of V2O3 in the corundum struc- 
ture. The coordination polyhedron of the vanadium atom is 
an octahedron of oxygen ligands trigonally distorted along 
the C"^ axis (i.e. two opposing planes of the octahedron are 
"squeezed" together). The point group of a vanadium site 
is Dsd. Moreover, inversion symmetry along the closest V-V 
bond is broken which leads to an onsite mixing of V 3d and 
V 4p states. 



lator transition upon doping or applying pressure^—. 
In the present work we not only provide all details of the 
previous calculation, but we also perform a completely 
new analysis for LD-XAS of V2O3. Last but not least, 
we exploit the XAS analysis of this important example of 
a metal insulator transition in a multi-orbital system in 
order to provide a general description of the K-pre-edge 
characteristics, which are by no means restricted to the 
case of V2O3. 

The paper is structured as follows: In section [III we 
recapitulate crystal and electronic structure of V2O3. In 
Section inil we give a pedagogical introduction to the the- 
oretical background for the absorption process with em- 
phasis given to the specifics of the K-pre-edge absorption. 
In the following section IIVI results of the density func- 
tional and configuration interaction are presented and 
compared to experimental data. 



II. CRYSTAL AND ELECTRONIC STRUCTURE 
OF V2O3 

At ambient conditions V2O3 is metallic and crystallizes 
in the corundum structure with four vanadium atoms in 
the primitive unit cell. The conventional and the primi- 
tive unit cells are sketched in Fig. [1] where it can be seen, 
that the vanadium atoms form "pairs" oriented along the 
crystallographic c-axis. Upon cooling below 150K, anti- 
ferromagnetic order sets in and the system becomes in- 
sulating; the onset of the magnetic order is accompanied 
by a transition from corundum to monoclinic structure. 
In addition, the system can be driven through (or further 
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FIG. 2. (color online) Level splitting of the vanadium t2g- 
states: due to the trigonal distortion (sketched for the aig 
case with the arrows as a compression along the threefold 
C^-axis of the octahedron) the t2^-states split up in an 
doublet and an aig singlet. From LDA results we estimate 
approximately 0.3 eV for this splitting. The plotted orbit als 
are spherical harmonic functions which display the symmetry 
of the states. 

away from) the Mott MIT by doping with chromium or ti- 
tanium or the application of external pressure. Above the 
Neel temperature the crystal structure does not change 
symmetry as a function of pressure or doping. How- 
ever, upon Cr doping a first order isostructural metal- 
to-insulator (PM-to-PI) transition takes place (termi- 
nating in a critical point arround ^ 400K) which evoked 
several theoretical attempts to describe this MIT in the 
spirit of a single band Hubbard model. While the MIT 
is associated to changes in the lattice structure and the 
atomic position o^^i^^ , it is important to notice that X- 
ray diffraction showed that for a given temperature the 
structure within one phase does not change upon dop- 
ing^ (but only with - for complementary infor- 
mation from photoemission see Ref. [25l ). It was later 
also observed by Park et at. with vanadium L-edge x- 
ray absorption spectroscopy that this holds also for the 
electronic ground state of the system (see Table 1 of Ref. 
[HI). Often it is assumed, that doping and pressure can 
be seen as equivalent routes through the phase diagram. 
Recently, however, experimental evidence and theoreti- 
cal analysis suggest to discard the common wisdom of 
the pressure/doping equivalence^. 

The electronic configuration of atomic vanadium is 
[Ar]3d^4s^ which means, that in the three- valent oxi- 
dation state one finds a 3d^ configuration realized. In 
the corundum type structure the vanadium atoms are 
coordinated by oxygen ligands in a trigonally distorted 
octahedral fashion (left hand side of Fig. [1]). Moreover, 
the crystal field breaks one significant point symmetry on 
the vanadium sites, namely inversion in the c-direction. 
This is related to the aforementioned V-pairs resulting 
in different distances of the neighboring vanadium atoms 
along the c-axis. As a result we find a slight onsite mix- 
ing of vanadium 3d- and 4p- states. While this effect is 



generally negligible for a discussion of the ground state 
composition, it will be of great importance for the se- 
lection rules of the polarization dependent XAS results 
which we discuss in the following sections. 

The major part of the crystal field is of cubic symmetry 
and, hence, splits the d-states into the lower lying i2g and 
the higher lying e^ states. The trigonal distortion acts 
like a compression along one of the three-fold axes of the 
octahedron (i.e. squeezing two opposite sides together). 
As a result the degeneracy of the lower lying t2g states is 
lifted and they are split into a single aig and an Cg dou- 
blet. In order to distinguish the two eg doublets, each 
of them gets an additional index indicating the type of 
bonds (either cj or tt) they form with the ligands: the high 
lying " cubic" eg (not split by the trigonal distortion) gets 
an additional " <j" and reads from now on whereas the 
lower one gets an additional tt and reads from now on ej. 
The level splitting, together with a plot of the respective 
angular part of the (atomic) wave function, is sketched 
in Fig. [2j Since the are higher in energy due to the 
crystal field, the two vanadium d-electrons populate the 
i2g levels. One of the crucial aspects concerning the un- 
derstanding of the MIT in the paramagnetic regime is 
the specific occupation of these i2g states. Experimen- 
tal evidence suggests^i^S^— that the ground state of the 
system should be described as a 5 = 1 state consisting 
of a mixture of aig and e^. Moreover, it is precisely the 
coefficients in the linear combination of aig and for 
the ground state which allow for a quantitative distinc- 
tion of the PM, PI, and AF phases. The XAS vanadium 
L-edge study of Park et al. explored the phase diagram 
by means of temperature and doping and summarized 
the respective ratios of aig and in table 1 of Ref. [l. 
As we mentioned earlier, their results turned out to be 
consistent with the X-ray diffraction data for the lattice 
of Robinsoni^, and showed that, within the PM- and 
within the PI - phase, there is no change in the ground 
state composition for different doping^^^. 



III. X-RAY ABSORPTION AND DICHROISM 
FOR THE VANADIUM K-EDGE 

In our recent work^^ we identified a robust ground 
state probe by calculating the ratio of the first two 
peaks in the K-pre-edge region. In Ref. [22| we discussed 
the origin of the ratio change in the case of a powder 
sample. However, a huge amount of extra information 
can be obtained by incorporating dichroic effects in our 
considerations. In fact, we will not only see that the 
dichroism of the K-pre-edge XAS fully supports the 
claim made for the powder spectra, but our theoret- 
ical results explicitely predict non-trivial interference 
between dipole and quadrupole transitions. The origin 
of such interference effects can be found in the earlier 
mentioned onsite mixing of vanadium d- and p-states: 
although the amount of p-character of the vanadium 
" d-orbitals" is small, the dipole allowed transition of the 
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Is core electron into this p-part makes the effect clearly 
visible in the K-pre-edge. However, the interference 
effect is not directly accessible in standard XAS since 
for these an average over all vanadium atoms in the 
unit cell is taken, leading to a symmetry conditioned 
cancellation. In fact, k-dependent absorption techniques 
could access the predicted intereference effects. 

Let us start by recalling the basic physics behind lin- 
ear dichroism in X-ray absorption. The reason for lin- 
ear dichroism and its relation to the occupancy of crys- 
tal field eigenstates can be understood intuitively even 
in a one-electron picture. Consider the following exam- 
ple: We want to promote a core 5-electron via a dipole 
transition to a valence state of p^-symmetr}^. Follow- 
ing Fermi's golden rule in the sudden approximation, 
the spectrum of such a process is proportional to the 
square of the integral {s\Dx^y^z\Pz) with dipole operator 
Dx,y,z (x^y^z)^ {s\ oc 1 and \pz) oc 2: in Cartesian coor- 
dinates. If the integrand is an odd function the integral 
vanishes. Hence, the only dipole transition operator for 
which the integral remains finite is the Dz cx z. The same 
holds, of course, also for the corresponding transitions to 
Px- or p^y-states, i.e., with x and y polarized light one 
can excite into the Px- or py-stdites respectively. Now, 
let the ground state of our hypothetical system be a sin- 
gle occupied p^-orbital and empty Px and Py orbitals. 
Hence, adding an electron to the p^-orbital will cost ad- 
ditional coulombic energy, so that the peak of the z- 
polarized spectrum (corresponding to the Dz transition 
operator) lies at higher energies compared to the Dx/y 
spectra showing also only half the intensity due to its 
occupancy. This means valuable information about the 
ground state of the system is encoded in the LD spectra. 
In other words: The polarization of the x-ray light adds 
more selection rules to the absorption process which can 
be used for a more detailed study of the ground state 
properties. 

Let us summerize now the most important transitions 
for the V i^-edge: The main edge consists mostly of 
dipole transitions of the s-core electron into the unbound 
V 4p-states, whereas transitions to the localized and cor- 
related 3d-states via higher order processes contribute to 
the pre-edge. We formally write the absorption process 
employing Fermi's golden rule as: 
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5{UJ — £^initial + £^final)- (1) 



This quantity is directly related to the spectral inten- 
sity. The transition operator describes the coupling of 
the photon to to our system and can generally be written 



^31.32. 



(2) 



where k and e are the light propagation and polarization 
vectors and r and p are the electron position and mo- 
mentum operators. Equation (|2]) can be understood in 



terms of different multipole transitions. Expanding e^^'^ 
in terms of multipoles gives: 



ikr 



J2 E ^'i2i + i)ji{kr)c^S{ek,<Pk)ciiH0r,<Pr), 



1^0 ra^-l 



(3) 

with ji being a spherical Bessel function of order /; /c, Ok^ 
and r, 6>r, (j)r the length and angular variables of the 
corresponding vectors k, and r, and finally Cm a renor- 
malized spherical harmonics Cm = \/47r/ (2/ + l)Ym'' . If 
one furthermore replaces the momentum operator p by 
the position operator r, using the commutation relations 
p = im/h[H,T] one finds for the transition operator: 

T.,k=Jo(fer)(r-e) (4) 
+ 3zji{kr){r • e)(r • k) + higher order terms 

(5) 

with r, k, and e denoting unit vectors in the direction of 
the position, momentum, and light polarization, respec- 
tively. 

In order to rewrite the full transition operator Te,k in 
terms of cubic-harmonics, we first define the multipole 
vector M consisting of the Dipole and Quadrupole oper- 
ators: 



M = {Dx,Dy,Dz,Qyz,Qxz,Qxy,Qx^-y^,Qz^) (6) 

The dipole operators can be expressed as: 
Dx = jo{kr)x 
Dy =jo{kr)y 

Dz=jo{kr)z. (7) 
The quadrupole operators read: 
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= ^6jl{kr)^/l/4{x'^ 




Qz^ 


= ^6jl{kr)^/3/4.{z^ 


-1/3) 



(8) 

In this basis we can define a generalized polarization 
vector taking into account only dipole and quadrupole 
terms, i.e., neglecting higher order terms, by means of an 
eight-component vector: 

^={^x, ^y, e^, (l/2)(e2^k2 + e^k^), (9) 

(l/2)( ^x)^z ~^ ^z^x ), (l/2)( ^x^y ^y^xj") 

^/l/4{exkx - eyky), y^iczkz - {l/3)e ■ k)). 

The full transition operator T^^k = M • e is then given by 
a dot product of ([6]) and (p!Q]) . Finally, the conductivity 
tensor can be computed as: 
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FIG. 3. (color online) Polarization dependent experimental data for the vanadium K-edge (left) compared to calculated LDA 
data (right). The color of the respective plot codes the geometry of the measurement, i.e., the transition operator, explained in 
the sketch on the left hand side. We can see the effects of the inversion symmetry breaking along the V-V bonds clearly in the 
insets, where we observe the splitting of spectra which should be equivalent by definition in a system without the symmetry 
breaking (see Eq. ()16p ). Note that the differences in the intensities of the largest peaks can be attributed to self absorption 
effects, which are not included in the calculation. 



Knowing the conductivity tensor, one can straight for- 
wardly calculate an absorption spectrum as —lm[e^-a-e]. 



Linear dichroism at the vanadium X— edge 

Let us now turn to the linear dichroism of the vana- 
dium K-edge. In Fig.[3]on the left hand side we show the 
experimental spectra (dots) of the V K-edge for three 
different polarizations of the x-ray light. In the upper 
inset we sketch the polarization of each spectrum with 
respect to the trigonal (local vanadium) reference frame 
(a, 6, c) (a and b form an angle of 120 degree, point to- 
wards the neighbouring V atoms in the plane and are 
both perpendicular to c which points in the direction of 
the out of plane V pairs). Since we would like to express 
the transition operators for each polarization in terms of 
the dipole and quadrupole operators ([7|) & (|8|) , we trans- 
late to Cartesian coordinates: 

a =x c =z (11) 

b=^y-\x y=-L(a + 26) (12) 

Now we can write down the dipole and quadrupole oper- 
ators for each orientation in terms of ([7|) and (|8|) 

Di, =D, Qb =Qyz (13) 

=\Dy + -D, Qg =g.2_^2 (14) 

/3 1 

=D, Qr =^Qxz + ^Qyz (15) 



Employing the crystal symmetry a ^ b we find straight 
forward 

Di, =D^ Qb =Qr (16) 

Evidently, while for dipole transitions we find the blue (b) 
and green (g) polarization to be equivalent, in quadrupole 
transitions blue (b) and red (r) polarization should be 
indistinguishable . 

However, when we inspect the experimental spectra 
closely we find these symmetries obviously violated: The 
lower inset of Fig. [3] (left hand side) shows a zoom of 
the pre-edge region where clearly all three spectra are 
different. Moreover, even in the main edge we find the 
"blue/green" dipole symmetry violated. As already men- 
tioned in section |TT] the reason for this violation can be 
understood intuitively with the sketch of the primitive 
unit cell (Fig. [T] right hand side). Along the connection 
line of two vanadium atoms in the primitive unit cell, 
i.e., parallel to the crystallographic c-axis the vanadium 
site lacks inversion symmetry, since the distances to the 
"upper" and "lower" V-neighbors is different. The con- 
sequence of breaking such a fundamental symmetry can 
be quite dramatic, as we will see. Namely, the V 3d- 
states are no longer pure eigenstates of the system and 
can mix with V 4p-states onsite. This means, that the 
formerly pure d-states get a tiny dipole moment and the 
quadrupole/dipole selection rule no longer holds strictly. 
This mixing will explain why the pre-edge is not entirely 
quadrupole in character and the main-edge not entirely 
dipole. Let us remark that this mixing is, in fact, small 
enough to be negligible for the formulation of the ground- 
state in terms of aig and occupation. In the absorp- 
tion experiment, however, the smallness of the mixing 
is compensated by the difference in the matrix-element 
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for dipole and quadrupole transitions. As a consequence, 
even a tiny amount of p-character is visible in the pre- 
edge, because it is dipole allowed: 



^|(^i|P|^f)| S{uj-AE) 



^ |(*5ri(«^l*4p) + bQ\^s,))f 5{co - AE) 



(17) 



where we defined AE = Ef — Ei as the energy difference 
of initial and final state. Now the spectral intensity of 
(fT7|) involves a dipole spectrum \{'^lf^\D\'^4p)f and a 
quadrupole spectrum 6^ |(^5s^®IQI^3c^)l^- However, due 
to the mixed character of the final state |^4p) + l^sd), 
both parts are entangled. A quantum interference term 
of the form 



ab{^sd\Q* 



(18) 



emerges which has to be taken into account explic- 
itly. In Fig. |4] we show a matrix representation of the 
spectra resolved in diagonal elements, i.e. pure dipole 
(rows/columns 1-3) and quadrupole (rows/columns 4-8) 
transitions and the off-diagonal interference terms. The 
different colors in each of the panels encodes the index of 
the four vanadium atoms in the unit cell. It is important 
to notice that when an average over all atoms is taken, 
as it is done in the present experiment, the interference 
terms cancel out exactly and are not directly accessible 
anymore. This is related to the fact that although inver- 
sion symmetry is broken for the pointgroup of a single 
vanadium site, it is not broken for the unit cell of the 
corrundum structure. Note that this would be different 
in the case of x-ray diffraction in which one measueres 
not the sum of —Im[e • ae] over all sites, but ^e^^'^^Fi. 
This means a phase dependent sum over the site depe- 
nent scattering tensor, whereby F = uja. In this case the 
inteference terms will not cancel out due to summation^^. 

The main edge, i.e., where excitonic features and the 
many-body coupling to the core hole are negligible, can 
be well described within the LDA. Since the ab initio 
LDA calculation incorporates the symmetry breaking in 
the structural input, the information about the resulting 
d-p mixing is already included in the eigenstates of the 
LDA Hamiltonian. 

In order to calculate the XAS spectrum we have to 
consider only the muffin tin sphere around the vanadium 
atom Rmt since the V Is core hole wave function is zero 
for r > RmT' We calculate projections of the LDA wave 
functions on the V 3d and V 4p subspace within the 
sphere: 



where we neglected the energy dependence of the ra- 
dial part of the wave functions, and the energy depen- 
dent Vi^rniuo) stems from the projection. The square 
of Vi^rn{uj) is in fact the so-called partial density of 
states. In terms of (p!9|) the integrals of (fT7|) separate into 
three parts: i) partial DOS '^sd{4p)(^) "interference 
terms" P3(i(cj) • V/^piuo) ii) Radial transition probabilities 
of the form 



/ r''dr'rRl^{r)R^p{r) 
Jo 



for the dipole part, and 

t-Rmt 



PltMT 

/ r^dr-r^Rl{r)Rsd{r) 
Jo 



(20) 



(21) 



for the quadrupole part. And finally iii) angular integrals 
like 



(22) 



(19) 



where O]? is the angular part of the specific transition 
operator. The angular integrals yield the selection rules 
for the transition. 



IV. LDA AND CI RESULTS 

In Fig. [3] on the right hand side, we present LDA spec- 
tra for the each transition symmetry. The overall agree- 
ment for the main edge is very satisfying. Especially 
effects of the 3d/4p mixing are captured well by the cal- 
culation and yield interesting features. As a signature of 
this mixing we find a splitting of the three spectra (see 
insets) in the experimental and calculated data, which 
would not be present in pure dipole/quadrupole transi- 
tions. 

In the pre-edge region the LDA calculation misses, ob- 
viously, the excitonic features (the first two peaks) which 
involve many-body states beyond LDA. Nonetheless, the 
LDA captures well the correct trends in the pre-edge 
region around 5470eV where transitions to the e^, the 
least correlated and most itinerant among 3d-states, are 
located. In order to use the "rule of thumb" from Ref . [22I. 
i.e., the ratio of the first two excitonic peaks (A and B) as 
probe of the a|e^, -f /3|e^, aig) (a^ + /^^ = 1) ground 
state, we have to include the 3d-4p mixing in our clus- 
ter calculation. This can be done by taking this mixing, 
i.e. the overlap integrals of V 3d and V4p, as a fitting 
parameter into account. The proper symmetry (but not 
the quantitative values!) of the potential we obtain from 
a Madelung calculation. As a technicality we remark 
that the t2g and spectra were calculated separately 
and summed up later. This was done in order to account 
for the big differences we expect for the broadening of the 
excitations due to the much stronger non-local hybridiza- 
tion of the orbit als with the oxygen ligands. Moreover, 
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FIG. 4. (color online) Conductivity tensor — Im[cr] (see Ea. (|lQp ) of V2O3 for dipole and quadrupole transitions at x-ray K-edge 
energies. Each panel shows a spectral function (for each of the four V atoms in the unit cell), representing an element of the 
conductivity tensor. The measured spectrum for a given experimental geometry is a linear combination of the tensor elements 
given by the generalized polarization (Eq. ^TO\i ) and the conductivity tensor as — Im[£* - a ■ e]. 



we calculated the spectra for all four vanadium atoms in 
the unit cell and averaged them. While locally the vana- 
dium atoms are equivalent, they differ with respect to an 
external frame of reference given by the propagation and 
polarization vector of the x-ray light. 

In Fig. [5] we report the results of our CI together with 
experimental data for the PM and PI phase. On the left 
hand side of the figure we show a "map" of simulated 
spectra for different values of a^. Overall we observe 
that we can capture the trend of the spectrum hierarchy 
also in the cluster calculation quite well. This time, how- 
ever, the itinerant states are beyond the basis set of our 
CI calculation of a small cluster and we cannot expect 
quantitative agreement for energies in the region of the 
Cg or higher. On the other hand, we now capture the 
excitonic peaks A and B correctly which we need as a 
ground state probe. Let us remark at this point, that 
the e^-states can be captured from both sides, LDA and 
CI, only qualitatively. In fact, they can be understood to 
represent a kind of " intermediate" regime - neither truly 
localized nor totally itinerant. The panel on the right 
hand side shows the color coding for the three polariza- 



tion directions. 

In the two middle panels of Fig. [5] we show the exper- 
imental spectra compared to the theoretical CI results 
which we obtained with the same mixing parameters we 
used for the isotropic spectra, i.e., very close to our pre- 
vious DMFT results^^. The agreement between experi- 
ment and CI calculations is quite satisfying. Moreover, 
we observe that the change of the ratio between peak A 
and B, i.e., our most important probe, going from PM to 
PI is much more pronounced in the selected LD spectra 
compared to the (averaged) powder spectra (see Ref f22l) : 
In the isotropic case the main spectral weight transfer 
responsible for the A/B ratio change was the one from 
peak B to higher energies, yet, in the LD data additional 
selection rules affect also peak A more directly, which 
has spectral weight almost proportional to the amount 
of Icg^aig) character in the ground state. This remark- 
able effect is, of course, due to LD selection rules and it 
works in a way just the same as the toy-example of the 
Pz ground state discussed at the beginning of Sec lIIIi It 
can be understood best in the limit of a pure quadrupole 
transition, where we neglect the mixing of the V 3d- 
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FIG. 5. (color online) Polarization dependent spectra in the 
pre-edge region of the vanadium K-edge. In the panel on the 
left hand side we show a "map" of theoretically calculated 
full multiplet spectra for different compositions of the ground 
state a\eg,eg) + /3|e^,ai^) (with + = 1). Employing 
this map we find the best agreement with the experimental 
data (middle panel) for values of — 0.5 for the PM phase 
(200K) and = 0.8 for the PI phase (300K). This provides 
a further confirmation that the LDA+DMFT densities^ are 
a good starting point for the CI calculation. 



states with the V 4p-states. In this limit the first peak is 
for our quadrupole transition operators, in fact, displayed 
exclusively in the |e^, ai^) part of the spectrum and com- 
pletely absent in the |e^, ep part. The transitions associ- 
ated with this peak is the \e^^aig) ^ e^, ai^, Iscore)- 
Note, if the ground state is |e^,aic,) the Is-core elec- 
tron has to be placed in the free level which can be 
done. If, however, the ground state is |e^,ep the Is- 
core electron has to be placed in the free aig level - but 
this is impossible with the specific transition operators 
of the displayed spectra, since all the associated matrix 
elements are integrals over odd functions and thus zero. 
Hence, we are directly sensitive to the aig occupation 
of the ground state. Let us stress that this argument 
holds strictly only in the abovementioned assumption of 



no 3d-4p mixing. In the real situation, and also our full 
calculation, this mixing weakens this effect slightly (the 
first peak is displayed also in the |e^, part of the spec- 
trum due to the d-p-hybridization) but nonetheless it is 
still clearly observable. In conclusion, the experimental 
LD data of Fig. [5] presents further evidence for the ro- 
bustness of our ground state probe and their theoretical 
simulation is consistent with the ground state parameters 
we found for the isotropic case. 



V. SUMMARY AND OUTLOOK 

By hands of the well known transition metal oxide 
V2O3 we have shown that XAS for the transition metal 
K-edge can provide valuable information on the elec- 
tronic structure of the d-orbitals, even though these are 
only accessible by quadrupole transitions and hence only 
form a preedge, which has an order of magnitude less 
weight than the main dipole- allowed oxygen edge. As a 
matter of course this low intensity is a disadvantage in 
comparison to the Vanadium L-edge. However, the hard 
x-rays of the K-edge can penetrate the diamond anvil 
cell so that XAS experiments are possible under pres- 
sure, e.g., to study the pressure induced Mott-Hubbard 
transition in Cr-doped V2O3. 

Our theoretical analysis has been based on LDA and CI 
calculations. The CI better accounts for local Coulomb 
interaction and excitonic effects, both of which are of im- 
portance for the d-states of the preedge, whereas LDA is 
more suitable for the more extended oxygen states. We 
elaborated the formahsm for the dipole and quadrupole 
transitions for different polarization directions of the 
XAS light, and elucidated on the quantum mechanical 
interferences between both. The dipole transitions are 
possible since there is a small, non-negligible 4p character 
of the predominantly 3d eigenstates in a transition metal 
oxide. Because of this there are two transitions paths 
of comparable magnitude from the core state to the pre- 
dominantly 3d eigenstates: (i) the quadrupole transition 
to the 3d part of the eigenstate and (ii) a dipole tran- 
sition to the mixed-in oxygen 4p part of the eigenstate. 
Because of quantum mechanics, we have (iii) interfer- 
ences between the two paths. To the best of our knowl- 
edge, these interference processes have not been studied 
in detail before. Let us note that for a lattice with inver- 
sion symmetry these interference terms cancel. However, 
for V2O3 there is locally no inversion symmetry for each 
Vanadium atom since there is a neighboring V atom in 
one direction along the c axis but not in the opposite di- 
rection. Hence there is a dipole-quadrupole interference 
term in the conductivity tensor (FigH]) for each Vana- 
dium atom. The whole crystal with 4 Vanadium atoms, 
on the other hand, is inversion symmetric; and accord- 
ingly the interference terms cancel in the average over all 
four vanadium spectra, i.e., in a XAS experiment. This is 
different for an x-ray diffraction experiment where, due 
to the k-dependence, the four atoms contribute with a 
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different phase. 

Of particular importance is the study of the polariza- 
tion dependence of the XAS spectrum, i.e., linear dichro- 
ism. This is necessary for the aforementioned interference 
effects, and also to extract information on the orbital oc- 
cupations, i.e., in the case of V2O3 on how many and 
aig states are occupied and how this changes at the pres- 
sure induced Mott-Hubbard transition^^. 

With our study, we demonstrated the usefulness of an- 



alyzing preedges in XAS with quadrupole transitions. We 
have highligted in detail the difficulties which are en- 
countered in such analysis and we showed how to solve 
them. Hence, we hope to trigger more XAS experiments 
of this kind as well as related electron energy loss spec- 
troscopy (EELS) experiments for which dichroism mea- 
surements recently became possible^^. We also propose 
x-ray diffraction experiments for V2O3 to experimentally 
prove the dipole-quadrupole interference. 
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